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THE  BEARING  CAPACITY  OF  A FOOTING  ON  A SOIL  (PLANS- STRAIN)1 

2 

by  R.  Tk  Shield 


Brown  University 


The  purpose  of  this  report  is  to  outline  a method  of 
showing  that  the  Prandtl  stress  solution  [l]^  to  the  plane- 
strain  problem  of  a flat  rigid  punch  (or  footing)  is  valid  for 
soils  with  cohesion  whose  angles  of  internal  friction  are  less 
than  75  degrees*  Since  a kinematically  admissible  velocity 
field  can  be  associated  with  the  Prandtl  stress  solution  [2], 
limit  analysis  [3,1*]  shows  that  the  Prandtl  value  [1] 

p = c cot  cp  [exp(x  tan  9)tan2(^5  +1L)  - 1]  (1) 

is  an  upper  bound  for  the  collapse  value  of  the  average  pres- 
sure over  the  punch,  where  c is  the  cohesion  and  9 is  the  angle 
of  internal  friction  of  the  soil*  A statically  admissible  ex- 
tension of  the  Prandtl  stress  field  into  the  rigid  region  is 
found  here,  so  that  the  value  (1)  is  also  a lower  bound  and 
therefore  the  true  value  of  the  average  pressure* 

In  Fig.  1,  OE  is  the  center-line  of  the  punch  OA 
indenting  the  surface  OD  of  the  soil.  The  region  OBCD  is  com- 
posed of  the  Prandtl  stress  field  of  two  regions  of  constant 
state,  OAB  and  ADC,  and  a region  of  radial  shear  ABC.  The  Prandtl 

■ . mi 

1The  results  presented  in  this  paper  were  obtained  in  the  course 
of  research  sponsored  by  the  Office  of  Naval  Research  under 
Contract  N7onr-3‘'801  with  Brown  University, 
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Research  Associate  in  Applied  Mathematics,  Brown  University, 

^Numbers  in  square  brackets  refer  to  the  bibliography  at  the  end 
of  the  paper. 
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field  is  continued  to  the  line  of  stress  discontinuity  BSP  which 
is  as  yet  unspecified.  The  stress  field  below  this  line  is  de- 
termined from  the  following  conditions,  (i)  dx  and  ay  are  func- 
tions of  y only  arvd  x only  respectively  and  is  2ero,  satisfy- 
ing the  equilibrium  equations  and  the  symmetry  condition  across 
OE.  (ii)  In  the  immediate  neighborhood  of  the  line  of  discon- 
tinuity, the  material  is  in  a plastic  state  of  stress.  These 
conditions,  together  with  the  jump  conditions  [5]  across  the  line 
of  discontinuity,  are  sufficient  to  determine  the  stresses  in 
the  region  BEFSB  and  to  determine  the  line  BSP.  Figure  2 shows 
the  Mohr's  circles  for  the  states  of  stress  on  either  side  of 
the  line  at  a typical  point  S.  It  is  found  that,  in  region 
BEFSB,  Ox  is  a monotonic  algebraically  decreasing  function  of  y 
and  dy  is  a raonotonic  algebraically  increasing  function  of  x. 

The  yield  condition  is  nowhere  violated  in  this  region  if  9 is 
loss  than  75  degrees. 

A perfectly  plastic  material  may  be  considered  as  a 
soil  for  which  9 — ^ 0.  Designating  the  shearing  stress  on  a 
plane  of  slip  by  k,  it  follows  from  the  above  that  the  value 
(2  + n )k  is  a lower  as  well  as  an  upper  bound  for  the  collapse 
value  of  the  average  pressure  in  the  indentation  of  such  a mater- 
ial. A somewhat  similar  method  has  been  used  by  Bishop  [6]  to 
show  the  validity  of  the  "complete"  solution  to  the  v-notched 
bar  problem  for  a perfectly  plastic  material. 
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APPENDIX 

The  soil  is  assumed  to  be  a plastic  material  In  which 
slip  or  yielding  occurs  in  plane  strain  when  the  stresses  satisfy 
the  Coulomb  formula  [7] 

2(°x  + °y)sln  9 + “ °y)2  + ^xy}^2  “ c cos  9 = °*  <2) 

This  equation  and  the  two  equations  of  stress  equilibrium  (in 

which  the  weight  of  the  soil  is  neglected)  form  a hyperbolic 

system  of  equations  for  the  determination  of  the  stresses  dx,  <jy 

and  x . The  two  characteristic  lines  are  inclined  at  an  angle 
xy 

n/4  + <p/2  to  the  direction  of  the  algebraically  greater  principal 
stress  and  will  be  called  the  first  and  second  failure  lines, 
with  the  convention  that  the  direction  of  the  first  failure  line 
is  obtained  from  the  direction  of  the  algebraically  greater 
principal  stress  by  a clockwise  rotation  of  amount  nA  + 9/2. 

The  anglo  of  inclination  of  the  first  failure  line  to  the  x-axis 
is  denoted  by  ©. 

It  is  convenient  to  put 
( o,  - d, ) 

p = — > 0,  (3) 

2 sin  9 

whore  0^  and  Og  (d^  < d 2 ' are  the  principal  stresses.  From  (2) 
and  (3)  it  can  be  shown  [8]  that 

dx  = - p [1  + sin  9 sin(2©  + 9) ] + c cot  9, 

d = - p [1  - sin  9 sin(2©  + 9)]  + c cot  9, 

y 

t = p sin  9 cos(2©  + 9). 
xy 
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The  equations  of  equilibrium  can  be  replaced  by  the  relations  [9] 


P cot  9 log  p + © - const,  along  a firs*  failure  line, 
cot  9 log  p - © = const,  along  a second  failure  line. 


(5) 


Across  a line  of  stress  discontinuity  in  a plastic 
stress  field,  the  normal  and  tangential  components  of  stress  must 
be  continuous  across  the  line  for  equilibrium.  The  equilibrium 
conditions  across  aline  of  discontinuity  separating  two  plastic 
stress  fields  £ and  £ can  be  written  [5] 


sin(©a  + ©^  - 2Q  + 9)  + sin  9 cos(©a  - 9^) 


(6) 


p cos(2©  - 2Q  + 9)  = p.  cos(20h  - 2Q  + 9),  | 
a a do  j 

where  Q is  the  inclination  to  the  x-axis  of  the  normal  to  the 
line  at  a current  point  of  the  line.  Subscripts  £ and  £ distin- 
guish the  values  which  p and  © assume  on  the  two  sides  of  the 
line. 


Referring  now  to  Fig.  1,  the  surface  AL  is  free  from 
traction  and  a = 0 in  the  constant  state  region  ACD.  Since 
© = x/4  - 9/2  in  ACD  it  follows  from  the  second  of  Eqs.  (4)  that 
p has  the  value  c cot  . 

p “ 

(1  - sin  9) 

in  the  region  ACD.  The  first  of  Fqs.  (5)  then  shows  that  at  a 
point  G on  the  first  failure  line  BCD, 


P = 


c cot  9 
(1  - sin  9) 


exp  1 2 tan  9 - | - a)  j 


1 N 
/ / 


where  a is  the  value  of  © at  the  point  G. 
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It  has  been  mentioned  above  that  the  material  just 
below  the  line  of  stress  discontinuity  BSP  is  assumed  to  be  in 
a plastic  state  of  stress.  We  shall  denote  by  £ and  the  two 
plastic  stress  fields  immediately  above  and  immediately  below 
the  point  S respectively,  where  S is  the  point  of  intersection 
of  the  second  failure  lino  AG  and  the  line  of  discontinuity. 
Since  the  failure  line  AGS  is  straight,  the  value  of  p,  pa,  in 
region  3 at  S is  also  given  by  Eq.  (7)  and  we  have  ©_  = a.  In 
region  J?  at  S,  ox  < so  that  ©^  = */4  - 9/2.  Also  the  normal 
to  the  1 _ne  of  discontinuity  at  S is  inclined  at  an  angle 
Q = n/2  - to  the  x-axis,  where  \J>  is  the  inclination  of  the 
lino  to  the  negative  x-axis.  Substitution  of  these  values  into 
the  first  of  the  jump  conditions  (6)  gives 


sin(^  + | + a + 29)  = sin  9 sin(^  + + a), 


and  the  relevant  root  of  this  equation  is 


1 3* 


a 

2 


4 

2 


9 

4 » 


where 


,7t  9 


sin  p = sin  9 sin(-^  + + a),  0 < p < 9. 


(8) 

(9) 


With  this  value  of  'i1,  the  second  of  tna  jump  conditions  (6)  gives 

C°s(jL  " a + P “ 


pb  = p 


a 


>•<5  i>  ’ 


or 

Pb 


c Jl+sin29-2  sin9Cos(^  + 2 +a-p)}  . 

= : —exp  I tan  9 -9  -2a)f  , (10) 

(1  - sin  9)sln  9 cos  9 2 J 


where  the  value  (7)  for  p.,  and  Eq.  (9)  have  been  used. 
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From  Eq,  (4),  the  non-zero  stresses  in  region  £ at  the 
point  S are  given  by 


n 


<3  ^ = - Pfc(l  + sin  9)  + c cot  9, 
<3  = - pfc(l  - sin  9)  + c cot  9* 

y 


> 

V. 


(11) 


For  equilibrium*  0 and  d are  taken  to  bo  functions  of  y only 
' a y 

and  x only  respectively  in  region  BEFSB.  The  values  of  c»x  and 

<Jy  just  below  the  line  BSF  are  known  from  Eqs,  (11)  and  (10)  30 

that  d and  d can  be  found  at  any  point  of  the  region  BEFSB. 
a y 

It  can  be  shown  that  p^  is  a monotonic  decreasing  function  of  a. 
It  follows  that  just  belcw  the  line  BSF,  d and  d are  increasing 
functions  of  a,  and  hence  that  in  region  BEFSB  d^  is  a monotonic 
decreasing  function  of  y and  dy  is  a monotonic  increasing  func- 
tion of  x. 


This  extension  of  the  Prandtl  stress  field  is  permiss- 
ible only  if  the  yield  condition  is  nowhere  violated  in  region 
BEFSB,  i.e.,  if  the  expression  on  the  left  hand  side  of  Eq,  (2) 
is  less  than  or  equal  to  zero  at  all  points  in  the  region.  Be- 
cause of  the  mono  tonic  character  of  the  stresses  dx  and  dy,  the 
yield  condition  will  not  be  violated  anywhere  in  the  region  if 
it  is  not  violated  at  the  point  E at  infinity  on  the  y-axis*  At 
the  point  E,  dx  has  an  (algebraic)  maximum  value  and  dy  has  an 
(algebraic)  minimum  value  and  these  values  are 


dx  = - 2 c tan(*  + ^), 

dy  = - c cot  9 jexp(n  tan  9)tan2(^  - 2)  - 1 J , 
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Substituting  these  values  into  the  expression  on  the  loft  hand 
side  of  Eq,  (2)  and  setting  the  resulting  expression  less  than 
or  equal  to  zero  gives,  after  reduction,  the  inequality 


exp(n  tan  9)  < 


tan^CjS 


!>■ 


The  inequality  is  satisfied  if  the  angle  9 is  less  than  an  angle 
which  lies  between  75  and  76  degrees.  Thus  the  yield  condition 
is  nowhere  violated  in  region  BEFSB  if  9 is  less  than  75  degrees. 

If  the  length  of  the  line  AS  in  Fig.  1 is  denoted  by 
r,  then  we  have 


or 


dr  = rd(a  + 9)tan(i|»  +0  + 9), 
dr  = rda  tan(^  + 


This  differential  equation  and  the  conuits  ori  r — AB  when 


a = - *A  - 9/2,  determines  the  line  jf  discontinuity.  As  a tends 
to  tiA  - 9/2,  the  line  tends  asymptotically  to  a straight  Line 
inclined  at  an  angle  nA  - 9/2  to  the  negative  x-axis. 
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Fig.  2 Mohr's  circles /for  states  of  stress  at  point  S 
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R-l  Attn:  R.  P.  Kroon,  Mgr.  of  Engineering,  ACT  Division 


8 


